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Abstract 

The relativistic Vlasov-Maxwell system of plasma physics is consid¬ 
ered with initial data on a past light cone. This characteristic initial 
value problem arises in a natural way as a mathematical framework to 
study the existence of solutions isolated from incoming radiation. Vari¬ 
ous consequences of the mass-energy conservation and of the absence of 
incoming radiation condition are first derived assuming the existence of 
global smooth solutions. In the spherically symmetric case, the existence 
of a unique classical solution in the future of the initial cone follows by 
arguments similar to the case of initial data at time t = 0. The total mass- 
energy of spherically symmetric solutions equals the (properly defined) 
mass-energy on backward and forward light cones. 


1 Introduction 

In a system of Cartesian coordinates the Vlasov-Maxwell 


system is given by 

dtf+p-V,f + {E+pxB)-Vpf = 0, (1.1) 

dtE-VxB = -j, dtB + V xE = 0, (1.2) 

V-E = p, V-B = 0 , (1.3) 

p{t,x)= / f{t,x,p)dp, j{t,x)= / pf{t,x,p)dp. (1.4) 


The Vlasov-Maxwell system models the dynamics of collisionless plasmas. We 
consider for simplicity a plasma consisting of a single species of particle. The 
unknowns are the particle density in phase-space, f = f{t,x,p), where is 

the momentum variable, and the mean electromagnetic field {E,B) = {E,B){t,x) 
generated by the particles. The expression 

^=—P— 
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denotes the relativistic velocity of a particle with momentum p. Units are chosen 
such that the mass and the charge of each particle and the speed of light are 
equal to unity. The simbol x denotes the usual vector product in . We refer 
to Pl El El m El ^1 ^1 for background on the Cauchy problem for 
the Vlasov-Maxwell system. Classical solutions of (n~Tl)-(n~ll satisfy the energy 
identity 

5te + V-p = 0, (1.5) 

where 

e{t,x)=J + + p{t,x) = Jpfdp + ExB. 

Integrating dni) one obtains the conservation of the total energy 

M{t) = j J + \p\'^ f dpdx + - j {\E\^dx = const. (1.6) 

Solutions of Vlasov-Maxwell also satisfy the continuity equation 

5tp + V-j = 0, (1.7) 

which upon integration leads to the conservation of the total (rest) mass 

N{t) = J J fdpdx = const. (1.8) 

The purpose of this paper is to set up a mathematical framework for the anal¬ 
ysis of solutions to the Vlasov-Maxwell system which satisfy the no-incoming 
radiation condition, that is 

rV2 r 

lim / / k-[E X B]{v—\x\,x)dSr{x)dv = 0, (1-9) 

J\x\=r 

for all Ui,U 2 SK, where dSr is the surface element on the sphere of radius r 
and k = x/r is the unit normal on this sphere. This corresponds to the physical 
condition that the electromagnetic held carries no energy to the past null inhnity 
of Minkowski space, see EIE]. 

Solutions of Vlasov-Maxwell isolated from incoming radiation were hrst stud¬ 
ied in El ■ The result of E] is that such solutions exist globally in time for small 
data of the Cauchy problem (i.e., data at time t = 0). However the Cauchy 
problem is not a natural framework to generate solutions isolated from incom¬ 
ing radiation. In fact, since the no-incoming radiation condition is imposed at 
t—>-—oo, there is no meaningful notion of local isolated solution with data at 
t = 0. Therefore, in the framework of the Cauchy problem, one can only prove 
the existence of global (or semiglobal) solutions which satisfy (ll.bll . This re¬ 
quires the use of uniform in time a priori estimates, which are not available in 
general for non-linear problems. 

A more natural setting for the study of isolated solutions is the initial value 
problem with data on a surface which cuts past null infinity. Examples of 
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such surfaces are past light cones and backward hyperboloids. This paper is 
concerned with the first case. 

Another motivation for studying the initial value problem with data on a 
past light cone comes from physical grounds. The initial data correspond to the 
outcome of an experimental measurement on the state of the physical system at 
the present time; the existence of a unique solution with the given data assures 
that the outcome of any future measurement is predicted by the theory. If 
this physical interpretation of the initial value problem is accepted, then it is 
clear that the initial data for relativistic models, such as the Vlasov-Maxwell 
system, should be given on a past light cone. In fact the set of events which are 
accessible to an observer at the proper time t = 0 lie on the past light cone with 
vertex on the world line of the observer at t = 0. The state of the system on the 
surface t = 0, on the other hand, cannot be measured, because these events form 
a spacelike hypersurface in Minkowski space. Such a discrepancy between the 
“physical” and the “mathematical” initial value problem has been sometimes 
discussed in the physical literature, see iziiHiiini and the references therein. 

In order to study the Vlasov-Maxwell system with initial data on a past light 
cone, we first rewrite the equations in the coordinates {v,x), where xSR.^ and 
u € R is the advanced time, which is defined by the condition that the surfaces 
V =constant correspond to the past light cones with vertex on the timelike curve 
|x|=0 (the world-line of the observer). Denote by f/\ = f/\{v,x,p),{E^,B^) = 
{E;^,B^){v,x) the particle density and the electromagnetic held expressed in 
these coordinates. They are related to the solutions of 1I.1II - (II.4II by f^{v,x,p) = 
f{v—\x\,x,p),{E/\,B/\){v,x) = {E,B){v—\x\,x) and therefore they satisfy the 
equations 

(1 -|-p- k)dvf/\ +P' Vx/a + {E/\ +px B/\) ■ Vpf/\ = 0, (1-10) 

dy{Er,-kxB^) = VxB/^-jf,, (l.II) 

dy{B^ k X E^) = — V X £//\, (1.12) 


dv{E/\-k)+ \7 ■ E/^ —P/\, ( 1 - 13 ) 

dyiB^-k)+V-B^ = 0, (1.14) 


where 


fc = l^, Pa{v,x)= J f/\dp, j,\{v,x)= J pf/\dp. 
Initial data are given at u = 0 and denoted by 

/A"(a^,p)=/A(0,x,p), A“(x) = £;a(0,x), B“(x) =Ba(0,x). 


(1.15) 


Later we shall discuss the equivalence of the system above with the evolution 
equations (n~Tni) - (n~Ta and a set of constraint equations on the initial data. 

In this paper we are interested in the question of existence and uniqueness 
of classical solutions in the future (i.e., for v€ [0,oo[) which match the initial 
data at v = 0. Obviously, one cannot expect (in general) that a unique solution 
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is determined by initial data at = 0, since the intersection between the initial 
surface and the domain of dependence of the solutions on a space-time point in 
the future is not a compact set. However it turns out that the Maxwell equations 
II. 1111 - 11.1411 have indeed at most one solution for given data at z; = 0 

provided the no-incoming radiation condition is satisfied. This suggests that 
the solutions we seek to the initial value problem with data on a past light cone 
should be restricted to the class of solutions isolated from incoming radiation. 

This paper is organized as follows. In Section|21we prove some general prop¬ 
erties of smooth solutions to the system 11.1011 - 11.1511 . The results of Section 
13 are conditional, as they assume the existence of global classical solutions. In 
Section 13 we also discuss the relation between the conservation laws satisfied by 
solutions with data on a past light cone and solutions with data at t = 0. Note 
in fact that for solutions with data on a past light cone, the conservation of 
the total mass and of the total energy are not obvious. In Section |3 we prove 
global existence and uniqueness of spherically symmetric solutions. This result 
is obtained by adapting to our case the proof of global existence for the Cauchy 
problem given in mini. In spherical symmetry the magnetic field vanishes 
identically (if decay at infinity is imposed) and the Maxwell equations reduce to 
the Poisson equation for the electric held. Hence there is neither incoming nor 
outgoing radiation in spherical symmetry. We will show that, as a consequence 
of the absence of radiation, spherically symmetric solutions satisfy the conser¬ 
vation laws iiTini . nn and that the total mass-energy equals the mass-energy 
on the past light cones and on the future light cones. 

In a subsequent publication the results of this paper will be extended to the 
Nordstrom-Vlasov system (see ^ for a derivation of this model). While it is easy 
to generalize the formal analysis of Section 2 below to the Nordstrom-Vlasov 
system, the proof of global existence and uniqueness of spherically symmet¬ 
ric solutions is different and considerably more involved since the Nordstrom 
scalar held equation remains hyperbolic—and so radiation propagates—also in 
spherical symmetry. 


2 The initial value problem with data on a past 
light cone 

An assumption on the initial data which will be made throughout is that 
0s^/“gC'^^(M3xK3), £;“,H“gC'2(IR3) 

and we dehne 

i?o = inf{i?: /™(x,p) =0, |a;| > R^pGM.^}. 

Hence = 0, for |a;| > i?o- In this section we study several properties of global 
solutions satisfying the regularity condition 

f^GC\[0,oo[xR^xR^), Aa,HaGC'1([0,oo[xIR3) 
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and so they are solutions of 11.1011 - 11.1511 in a classical sense. We also assume 
that /a has bounded support in the momentum, precisely 

7^a('c) = sup{|p|: 7 ^ 0, 0 < s ^ u, < oo, VuGM. 

In particular, all the integrals in the momentum variable in the sequel are under¬ 
stood to be extended over a compact set. We split the analysis in two different 
subsections. 


2.1 The Vlasov equation 

We start by pointing out some basic properties of /a. Note the estimate 


1+p-k = 1-1- 


P-k \p\ 

\/i + bP ^ V^ + \p\ 


1/2 . 

xAWCvTTh^+H) " i + bP’ 




( 2 . 1 ) 


hence when the support in p of /a is bounded, the equation iurnii is equivalent 
to 

Qv/a H-H-( V^ld" IpP E/\ +p X Ba') • Vp/a = 0, (2-2) 

Po Po^ ) 

where po is defined by _ 

po = ^/l+]p^+p-k>0. 

The characteristics of the differential operator in the left hand side of (E3 are 
the solutions of 


x = —, p= — (\/1 -|-|pPBa-I-pxBa) 
Po Po ^ ^ 


(2.3) 


and we denote by {X,P){s,v,x,p), or simply {X,P){s), the characteristic sat¬ 
isfying {X,P){v) = {x,p). Since the particle density /a is constant along these 
curves, we obtain the following representation formula for the solution of the 
Vlasov equation: 

f^iv,x,p) = f^{{X,P){0,v,x,p)). (2.4) 

In particular /a remains non-negative for all times and ||/a(u)||oo ^ ||/a”IIoo- In 
the next lemma we estimate the a:—support of /a. 


Lemma 1 For all v^O, 


f^{v,x,p) 


0, for\x\^Ro + -v. 


Proof: For all 0 ^ s ^ u we have, by the first equation in (IO|l . 


|a:| = |X(0)| + 


P{t)-K{t) 


'o v'l + |P(T)|2+P(r).iF(r) 


dr, 
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where K = X/\X\. Let [0,t)]=X UX+, where 

X- = {r G [0, z;]: (P • X) (r) 0}, X+ = {r G [0, z;]: (P • X) (r) > 0}. 

Thus, using y^l + |pp >p-k, 


1^1 < l^(o)|- 


PK 


li+ y/T+\P\^ + P-K 


dr 


< |X(0)| + imeas(X+)<|X(0)| + iw. 


Since |-^(0)| ^Rq in the support of /a, the lemma is proved. □ 

We shall now derive the conservation laws satisfied by the solutions of ES)- 
A straightforward computation reveals that the right hand side of the system 
(toil . i.e., the vector 


F{v,x,p) = 


JL J_ 

Po'po 


^\/l + IpP E/\ +P X Bf^ 


satisfies 


[V(,,p)-P](s,A(s),P(s)) = -^log(l + P(s)-iG(s)), (2.5) 

where P = P/y/l + P^. In fact, each side of 12.511 equals, along characteristics, 

I ■ (fc - (F■ fc)p) - (P X fc) • Pa) 

\x\ v^l+|p|2 V J 


(l+p-kY 
From 12.511 we deduce 


det 


d{X,P){s) - 
d{x,p) _ 


1+p-k 

l + P(s)-iG(s)' 


Hence using (Ell) the next lemma follows. 

Lemma 2 For any measurable function Q : K —> K., 


J J Q{fA){^+P' k)dpdx = const. 


In particular, by choosing Q{z) =z‘^, 1, 

\\{l+p-ky/‘>f^(v)\\L 9(183 xR3) = const. (2.6) 

The case g = 1 in 112.611 corresponds to the conservation of the (rest) mass on 
the past light cones. To be more precise, observe that PajJa defined in HI. 1511 
satisfy the equation 

dviPA+jA-k) = -X -jA- (2.7) 
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The latter can be proved either by using (trrni) or by a simple change of variable 
in We define the mass N/k^{v) on the past light cone at time v as 

N^(v)= limn/-,(v,r), n/^{v,r)= {p/^+j^-k)dx. 

Note that the function nA(w,-) is non-decreasing and so the above limits exists. 
By (I2.till . Nfy{v) = Nfy{Q), for all v^Q. The total mass of a solution, given by 
can be rewritten as 


N{t)= lim n(t,r), n(t,r)= / p{t,x)dx= / p/^{t+\x\,x)dx, t^O. 

./|a:|5Cr J\x\^r 

In the next lemma we prove a formula which relates the mass functions N(y) 
and N/\{v). 

Lemma 3 For all v^O, 

pv-\-r p 

n{v,r)=n^{v,r)- / ,x)dSr{x)dv'. 

JV J |a:|=r 

Proof: Integrating between v and u-l- |a:| we get 


(Pa +Ja • k){v + \x\,x) - (pa + Ja • k){v,x) 
Integrating in the region |a;| we get 


>+\x\ 


'^■j/x{v\x)dv'. 


/ {PA+jA-k){v + \x\,x)dx= {p^+jA-k){v,x)dx 

J |ai|^r J lail^r 


Now we use the identity 
u+l^l 


|ai|^r 


( 2 . 8 ) 


u-|-|x| 


| 3 ;| J V 


^ ■ JaW ,x)dv' dx. 


V- 


r'V+\x\ r■v+\x\ 

/ JaW, x)dv'= {j^-k){vF\x\,x)+ / V ■ j^{v',x)dv'. 

J V J V 


Substituting into llOll and using the Gauss theorem proves the lemma. □ 
By Lemmanall characteristics of the Vlasov equation must cross the surfaces 
t = v + \x\=const. for all t^O in compact sets of x. This means in particular 
that no mass can be lost at spacelike infinity, which explains why the following 
lemma holds true. 


Lemma 4 For all v'^0, N{v) = IVa (v) = N/^ (0). 
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Proof: By Lemma ^ and Lemma|21we have n{v,r) =n/\{v,r), for r>2Ro + v. 
The claim follows by letting r—> 00 . □ 

Next we define the mass on the future light cone at time v as 

N'^{v)= lim n^(u,r), n^(u,r)= / {p^ — j'^ ■k)dx, 

dtxKr 

where 


p^{v,x)=pf,{v + 2\x\,x)=p{v+\x\,x), 
j'^{v,x)=j^{v + 2\x\,x)=j{v + \x\,x). 

By a change of variable in (lOl (o r in IE3) we have 

d4p'^-f-k)=-V-j'^, (2.9) 

Integrating (1^ in time between v and r; + 2|a;| and proceeding as in the proof 
of Lemma 13 we obtain 

pv-\-2r p 

n'^(z),r) = n/\(u,r) — / / j/\-k{v',x)dSr{x)dv'. (2-10) 

J V J |rc|=r 

Moreover, by (E3, satisfies the equations 

i9„n^ = — / j^-kdSr{x), — dr-n^ = — f p'^dSr{x). (2-11) 

J\x\—r J\x\—r 

The evolution of the mass on the future light cones is studied in the following 
lemma. 

Lemma 5 The function (v) is non-increasing, that is, 

{v 2 )^N^{vi), \lvi^V 2 - 


Moreover 

pV2 

{i) N'^{V 2 ) = {vi) iff lim / 

r^oo / 

J Vi 

pv-\-2r p 

(a) N'^{v) = N^{v) iff lim / / j/^-k dSr{x)dv'= 0. 

J\x\=r 

Proof: Integrating the second equation in itmii along characteristics we have 
{v 2 ,r — V 2 ) ivi,r — vi), for all t >2 ^ wi and r>V 2 - In the limit r^oo this 
implies that N'^ is non-increasing. The claim (i) follows by integrating in time 
the first equation in lITTTll on the interval [r'i,U 2 ] and letting r—^oo, while (ii) 
follows by lITTnil . again in the limit r—>oo. □ 

We remark that for solutions with data on a past light cone it is not obvious 
that N'^{v) is bounded. Moreover, even if bounded, it needs not to be constant. 


i 


\x\—r 


j'^ ■ k dSr{x)dv = 0; 
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If N'^{v 2 )<N'^ivi), for 0<i>i<z;2, the difference (vi) — N'^{V 2 ) measures 
the mass lost at future null infinity in the interval [^ 1 ,^ 2 ] of the advanced time. 
Finally, even if N'^ is bounded and constant it is not obvious that it must equal 
since the limit condition in (ii) of Lemma |31 might not be satisfied. 

In the next lemma we show that a sufficient condition for the limits in (i) 
and (ii) of LemmaElto be zero is that the momentum support of /a is bounded 
uniformly in u G K., as this condition implies that no particles can reach future 
null infinity. 

Lemma 6 Assume V/\{v) ^D, for all v^O and for some positive eonstant D. 
Then, for all vi,V2,v^0, 

N'^{v2) = N^{vi), N'^{v) = N^iv). 


In particular, by Lemma^^ N'^(v) = N(v) = N/\(v) = N/\{0), for allv^O. 
Proof: By the assumption, 


,- Vl + D^, , Vl + L>2 

\/l+bP5s-- \P\^ -- (p-k) 

in the support of /a and so, as in the proof of Lemma ^ 

P{t)-K{t) , _ D 


N<|X(0)|- 


/i+ a/ 1 + |P(r)p+P(r) •iL(T' 


-dT<Po + 




for all (x,p) G supp/a(u), where = {tG [0,u] : (P-iL)(T) > 0}. This implies 
that f^{v,x,p) = 0, for |a;|^Po + au and aG[0,i[. Hence 

rv 2 r 

/ / j'^ ■kdSr{x)dv = 0, for r>{l — 2a)~^{RQ + av 2 ), 

Hi J\x\—r 
pv-\-2r p 

/ / j/s,-kdSr{x)dv' = 0 , ioT r> {l — 2a)~^(RQ + av). 

JV J \x\—r 


Lemma 4 concludes the proof. 


□ 


2.2 The Maxwell equations 

We now pass to study some general properties of the electromagnetic field 
(Pa,Pa)- First we show the equivalence of the Vlasov-Maxwell system with 
a set of evolution equations and a set of constraint equations on the initial data. 
An important point is that, in the present situation, there are more constraint 
equations than in the case of the Cauchy problem, since the initial data are 
given on a characteristic surface. Computing the vector product of inu with 
the unit vector k, subtracting and then using itmii we obtain 

fc X (V X Pa) — fc(V • Pa) + V x Pa — fc x Ja = 0. (2-12) 
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Moreover, computing the vector product of with k, adding om and 

then using we obtain 


V X B/-, + fc(V • i?A) — fc X (V X i?A) — Phk — jh= 0. (2-13) 

On the other hand, the equation itrra follows from cm and cm, whereas 
(tmil follows from and cm. Hence the whole set of the Maxwell 

equations is equivalent to the system composed by n~TTll - (n~llll and HTWi - 
(tmil . Clearly, since cm-Em are valid for all times, then they must be 
imposed at v = 0 in order to obtain a solution of the initial value problem, 
i.e., cm-iEm are constraint equations on the initial data. However these 
constraint equations are not totally independent. Let Wi, W 2 denote the left 
hand side of and cm, respectively. It is easy to verify the following 

identities: 

Wi = kxW2 + k(k-VxE^-V-B^), 

W 2 = —k X Wi + k{k ■ V X B/\ + V • E/\ — P/\— j/\ ■ k ). 

From this it follows that cm-(Em are equivalent to the equations 

V ■ B/\ = A: ■ V X E /\, A: • V x B/\ F V ■ E/\ = p/\ F Ja ■ k^ (2.14) 

together with one of the equations k x Wi = 0, or A: x W 2 = 0, that is 

A:x [(A:x VxBa)FVxF;a-A:xja]=0 (A:xIFi=0), (2.15) 

or 

A:x [VxHa-A:x VxF;a-Ja] = 0 (A:xtF 2 = 0). (2.16) 

The following proposition concludes our discussion on the reduction of the 
Vlasov-Maxwell system to a set of evolution equations and a set of constraint 
equations on the initial data. 

Proposition 1 The following assertions are equivalent: 

(V (/aj^'Aj-Sa) is a solution to the initial value problem for 

(2) {f^,E/^,B/^) is a solution to the initial value problem for 
and the initial data satisfy 

(3) {f^,E/^,B/^) is a solution to the initial value problem for \l.l(A) - kl.BA) 
and the initial data satisfy 

(4) (fA,E/\,BA.) is a solution to the initial value problem for J^1.1(A) - \1.1‘A) 

and the initial data satisfy and 

Proof: We already proved that (1)<;4>(2)^(3) and (3)<;4>(4). Thus it is sufficient 
to establish (4)=^(1). It is a simple exercise of vector algebra to show that 
(imil are satisfied for all times provided they are satisfied at time ?; = 0 and 
£^a,.Ba,Pa,Ja satisfy lll.llll . 11.1211 and C3. The latter holds in virtue of the 
Vlasov equation cm- Moreover, (trrni follows from cm and the second 
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equation in while itmii follows from srwi and the first equation in 

12.1411 . Thus is a solution of ll.10ll - ll.14ll and since iirnii - iiniii 

are satisfied at v = 0, then it is also a solution of the initial value problem. □ 
The no-incoming radiation condition in the coordinates (v,x) reads as in the 
following 

Definition 1 A global solution of is said to satisfy the no¬ 

incoming radiation condition (NIRC) if, for all vi,V2^0, 

l‘V 2 r 

lim / / k-[E/,,x B/,^]{v,x)dSr{x)dv = 0. 

J\x\=r 

A local solution in the interval V>0, satisfies NIRC if the above limit is 

zero for all Vi,V 2 G [0,14[. 

We impose NIRC only in the future, since our purpose is to study the initial value 
problem forward in time. Likewise we may introduce the concept of outgoing 
radiation as in |^. 

Definition 2 The outgoing radiation emitted by a (global) solution 

of the Vlasov-Maxwell system in the interval of the advanced time is 

given by 

r'^2 p 

'Sout(i’i,i’ 2 ) = lim / / k-[E^x B/\]{v-\-2r,x)dSr{x)dv, 

J\x\=r 

provided the limit exists. 

The energy identity in the coordinates {v,x) reads 

9„(eA+pA-fc) = -V-pA, (2.17) 

where 

eA = y \/i+1pF/a 2 l-^A 1^ + 2 l-®A l^j pA= J pfxdp-I E/\ X B/\. 

The identity ETtIi can be proved either by a direct calculation using the equa¬ 
tions (trrmi - itrra . or by a simple change of variables in dni). Next define 


mA.{v,r)= {e/^-\-p/^■k)dx. 

J\x\^r 

By 12.17I) . ttia satisfies the equations 

dvm/\ = — p^-kdSrix), 9«mA-l-9rmA= / ef^dSr{x). 

J \x\—r J |af |=r 

We define the energy Ma(u) on the past light cone at time v as 

Ma(u)= lim mA(t',’'). 


(2.18) 


The function mA(f,’) is non-decreasing and so the above limit exists. 
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Lemma 7 is a non-decreasing function: 


Ma(wi) yVi^V 2 . 


Moreover if the NIRC is satisfied then Ma(i') is constant for all v^O. 


Proof: For all vi ^ V 2 and r>V 2 — vi we have, integrating the second equation 
in 12.1811 . mA(r'2, r) m/\{vi ,rvi — V 2 ) and letting r—>oo we prove that Ma 
is non-decreasing. To show that M/\ is constant in the absence of incoming 
radiation, we use that, by the first equation in 


m^{v 2 ,r)-m/\ (wi,r) 




Jp- kf/\dpdSr{x)dv 



k- [i?A X B/s]dSr{x)dv. 


By Lemma n the first term in the right hand side vanishes for r>2i?o+'L'2, 
while the second term tends to zero in the limit r^oo by the NIRC. □ 

From Lemma 0 we obtain the following uniqueness theorem for the Maxwell 
equations. 


Lemma 8 = 0 is the unique solution of the homogeneous system 


3„(i5A — fc X Ra) — V X Ra =0, dv{B^-\-kx X E/^ = 0, (2-19) 

which satisfies the NIRC and the initial condition {E;^,B^){0,x) = 0. 

Proof: By Lemma 0 we have 

0 = 2[\E^\‘^ + \B^\‘^ + 2iE^x B^)-k] 

= \E^ ■ fcp -|- |i?A ■ + \Ea — kx i?Ap + |Sa -l-fc X £/Ap- 


Hence the solution is a plane wave propagating along the fc—direction, i.e., the 
vectors {E/^,B/^,k) form an orthogonal triad. It follows by 12.1911 that V x E/^ = 
V X i?A = 0 and so, by 12.1411 . V • E/^ — V • Ha = 0- The claim follows. □ 

By a standard interpolation argument we obtain 

Lemma 9 If the initial data are chosen such that M/\(0) is bounded and the 
solution satisfies NIRC, then 

II (Pa + Ja • ^){'(')IIl4/3(b3) ^ C'.^a(O)) Vw 0, 

where C is a positive constant which depends only ||/a”||oo- 
Proof: We write 


PA+jA-k= / (l-bp-fc)/Adp+ / (l+p-fc)/, 

J\p\^R J\v\>R 


dp 


p\>R 


< YII/“IIooH^ + H ^ J Pofxdp^C (^J po/a^ 

^ C{eA-IpA-kf^'^, 
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where in the second line we choose 


R = 


(^II/a°IIoo^ J Poffyd^ 


The claim follows. □ 

We shall now briefly discuss the relation between the the total energy and 
the energy on the past light cones. The total energy Oil can be rewritten as 


M(t) = lim m(t,r), m(t,r)= / e{t,x)dx= / e^(t+\x\,x)dx. 

J\x\^r J\x\^r 


By (ESIl, m('u,r) satisfies the equations 


9,,m = — 


/ |x|=r 


P • kdSr(x), 


( 2 . 20 ) 


9„m±9rtn = 


j {±e —p ■ k)dSr{x). (2.21) 

J |x|=r 


The right hand side of (OTll is non-negative in the -I- sign case and non-positive 
in the — sign case. 

Lemma 10 The total energy is constant, i.e., 


M{v2)=M{vi), Vwi,U2^0. 

Proof: Integrating lOTTl with the plus sign along the characteristics of dv + dr 
we obtain m{v 2 ,V 2 +r)^ m(wi,ui -|-r) for all vi^V 2 , which implies, in the limit 
r—>c», M{v 2 )^ M{vi). On the other hand, integrating 112.21 |l with the minus 
sign along the characteristics of dy — dr gives m{v 2 ,r — V 2 ) ^m(wi,r —ui), for all 
V 2 ^v\,r>V 2 and so, letting r^oo, M{v 2 )^M{vi). The claim follows. □ 
We emphasize that for solutions with data on a past light cone it is not 
obvious that M is bounded. If it is bounded, then, by Lemma^J it is conserved. 
In the latter case, however, the total energy and the energy on the past light 
cones might be different. To see this consider the equation 

/'V-\-r p 

m{v,r)=m^{v,r)— / p^-k{v',x)dSr{x)dv'. (2.22) 

J V ^ \x\—r 

The latter is obtained by integrating ETtIi in time from v to u-|- |x| and pro¬ 
ceeding as in the proof of Lemma 13 By Lemma Hand 112. 2 211 we have 

m(v,r) =mA(u,r)— / f k-[Eff*'x ,x)dSrix)dv', (2.23) 

J V J |ai|=r 

for r>2i?o + u, where is the field outside the support of the matter. 

Hence the answer to the question whether or not = M depends on the decay 
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of the solutions of iITTnii as r —> oo. As we shall discuss in Section^ the equality 
Ma = M holds for spherically symmetric solutions, as in this case the magnetic 
field vanishes identically. An interesting open question is whether Ma = M holds 
in general in the absence of incoming radiation. 

To conclude this section we study the evolution of the energy on the future 
light cones. Let 

e'^(v,x) = e/^,(v + 2r,x) =e(v + r,x), 
p'^(v,x) =pA,(v + 2r,x)=p(v + r,x), 
which satisfy the equation 

dv (e^ — • fc) = — V • p^. 

Now let 


m'^(r;,r)= / (e'^ —p'^ ■ k)dx, M'^{v)= lim m^(z;,r) 

J\x\^r 

and note the equations 

dym'^ = — f p'^-kdSr{x), dyvn'^ — drvn'^ = — f e'^dSr{x), 
J\x\—r J\x\—r 


(2.24) 


pv-\-2r p 

{v,r) =m/^,{v,r) — / p/\-k{v',x)dSr{x)dv'. (2.25) 

J V J \x\—r 


m 


' |a:|=r 


By (12.2411 - (12.2511 and the same argument as in the proof of Lemma|5lwe obtain 
Lemma 11 For all vi ^V 2 , 

Moreover 


pV2 r 

(i) M'^(v 2 ) = M'^{vi) iff lim / / p^ ■kdSr{x)dv = 0\ 

J\x\=r 

pV-\-2r p 

{ii) {v) = M/\{v) iff lim / / p/\ - k dSr{x)dv' = 0. 

J\x\=r 


’ |tc|=r 

The remarks on N'^ following the proof of Lemma [S] apply to as well. In 
particular, the difference M'^(vi) — M'^{v 2 ), when it does not vanish, measures 
the energy dissipated by the system to future null infinity in the interval [^ 1 ,^ 2 ] 
of the advanced time. By (ii) of T/emma TTTl this is the sum of two contributions: 
an energy lost in form of outgoing radiation by the electromagnetic field (as given 
in Definition mi and a kinetic energy carried by the particles, which is given by 
the limit 


lim 

r—>-oo 


c|=r 


P‘ k f/\{v + 2r,x,p)dpdSr{x)dv. 


As in the proof of LemmalHl the latter term vanishes if the momentum support 
of /a is uniformly bounded in r; S R, as in this case no particles can move to 
future null infinity. Given this interpretation, it is natural to identify M'^ as 
the analogue of the Bondi mass in General Relativity, see [2] . 
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3 Spherically symmetric solutions 

In spherical symmetry we have V x = V x B/\ = 0 and so, by the first equation 
in I2.14II . V ■B/\ =0. Under the additional boundary condition liiUr^oo^A = 0, 
this implies that the magnetic field vanishes identically. Moreover, by the second 
equation in iimi . 

k r 

Ea{v,x) = ^ / {pA+jA-k){v,r'y‘^dr' 

Jo 

= ^/ (Pa +Ja • k){v,y)dy, (3.1) 

the second equality being valid in spherical symmetry. By abuse of notation 
we use the same symbol to denote a spherically symmetric function in spherical 
and Cartesian coordinates. The Vlasov equation reduces to 

d^A + - ■ ^ Ea ■ Vp/A = 0. (3.2) 

Po Po 

In spherical symmetry the particle density is invariant under proper rotations in 
phase-space. This allows one to write /a = fA{v,r,w,q), where w = (p- fc) € R and 
q= |xApp ^ 0, see m However the Vlasov equation is more conveniently stud¬ 
ied in the coordinates {x,p). Note also the conservation of angular momentum: 
along characteristics, 

^Jxxp\^ = 0. (3.3) 

In the spherically symmetric case we have the following global existence theorem. 


Theorem 1 Let xR^) be spherically symmetric and satisfy 

F = inf{ |a; x pp : {x,p) G supp/™} > 0; (3.4) 

there exists a unique, spherically symmetric G C'^([0,oo[xR^ x R^) solution 
of ld.l\} - ld. iH) such that //^(0,x,p) =/™(x,p). Moreover, there exists a constant 
C >0, depending only on bounds on the initial datum, such that 


Va{v)^C. 


(3.5) 


Before giving the proof of Theorem ^ let us observe the following 
Corollary 1 For the solution of Theorem^ 

N'^{v) = N{v) = Na{v) = Na{Q), 

M'^(v) = M{v) = Ma{v) = Ma{Q). 
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Proof: The equality of the mass functions follows from Lemma El Since spher¬ 
ically symmetric solutions are isolated from incoming radiation, then M/\{v) is 
constant by Lemma [T) Setting = 0 and letting r^oo in (I2.2dll . we have 
M{v) = M/\{v). Hence it remains to show that M'^(v) = M/^{v) , for all v^O. 
For this purpose we use 12.2511 with B/y = 0, that is 

|•v+2r /• p 

m^(w,r) =mA(v,r)- / / / p-kf/^{v',x,p)dpdSr{x)dv'] 

Jv J\x\=rJ\p\^C 

as in the proof of Lemma El the integral in the right hand side of this identity 
vanishes for r large enough and letting r—^oo concludes the proof. □ 

The proof of Theorem ^ is formally identical to the proof of global existence 
for the Cauchy problem with data at time t = 0 given in Theorem II] (see ^3] 
for the case of two different species of particle). We shall sketch it for the sake 
of completeness, restricting ourselves to derive the main estimates which lead 
to the proof. Note however that the assumption is not made in ini Here 
the condition ea is used to ensure that the characteristics are in all the 
parameters. In fact, due to the presence of the unit vector k, the coefficients of 
the Vlasov equation are in general discontinuous at a: = 0. But thanks to 
and the conservation of angular momentum, the solution is supported away 
from the axis jxj = 0 and so it can be defined in a classical sense in terms of the 
characteristics. The assumption ei can probably be removed by passing to a 
weaker solution concept, but we shall not pursue this here. 

For the proof of Theorem ^ we denote by C any positive constant which 
depends only on the initial datum. Moreover we define 

Rmiu{v) = inf{|x|: f^{s,x,p) 7^0, 0 < s^ v.pGM.^}. 


By the conservation of angular momentum and (lOl . 




Vf 

'Pr.iv)' 


(3.6) 


Hence a bound on the momentum support of /a implies that the particle density 
vanishes in a neighbourhood of the axis jxj =0. Now observe that, by (ESI) for 


9=1, 


\E^{t,x )\^ 


Moreover, the bound (pa + Ja ’ k) ^ CV/\{v)^, LemmaEland Holder’s inequality 
imply 




2/3 


C 

r 


11 Pa + J A • fc 11 iCa ^A (w) < C'T’a (?^) . 
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Next define 


G{v,r) = — min 
J r 

It follows that G{v,-) & is increasing and \E/\{v,x)\ ^drG{v,r), for all v^O. 
Moreover, and since 'Pa(-) is non-decreasing, drG{vi,r) ^drG{v 2 ,r), for vi <V 2 - 
Splitting the integral at R= y/W^{CV one obtains 

_ / , \ 1/2 

G{v,Q)=-2^^{cV;,{vf/^) 

and therefore, for all ri,r 2 ^ 0, 

|G(r;,ri)-G(t>,r 2 )| < |G(t>,0)| < CV^{vf/\ 

Next we claim that 


iVA 

A2 


CV;,{vf/Ad\ v,r^0. 


(tt) there exists at most one UoS[0,oo[ such that ^|X(s)|=0 and if such vq 
exists then |^(s)| has an absolute minimum at s = vq. 

This follows because, along characteristics, 

ds Po ds po \ \x\ ) 

Since p-k is increasing and ^|N'(s)| > 0 (resp. < 0) for p• fc > 0 (resp. < 0), the 
claim (tt) is proved. Now observe that, along characteristics, 

ds Po 

Hence, for all u ^ 0 and sG [0,u] we have, denoting J+ = {tG [0,s]: ^ |X(r)| ^ 0}, 

Vl + P(s)2-Vl + P(0)2= r\E^{r,Xir))\^^^l^^dr 

Jo ^o\t) 

= j^^\E,,{T,Xlj))\^\X{T)\dr 
= j^^j^[G{s,\X{r)\)]dr. 

By virtue of (tt), either J+ = [si,s], for some 0 < si < s, or J+ = [0,s], or X+ is 
empty. In the first case we obtain 

x/l + P(s)2 < Vl + m' + G(s,|X(s)|)-G(s,|X(si)|) 

< V^l + P(0)2 + GiPA(u)®/® 
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and since this is true for all then + which im¬ 

plies Va. (v) < C. The other two cases lead to the same inequality. The bound 
on the momentum support of /a implies, by Id.hll . that the particle density is 
supported away from the axis |a;|=0. This allows one to define /a in terms 
of the characteristics and derive L°° estimates for its derivatives. A standard 
iteration scheme completes the proof of the theorem. 
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